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Recently, Gohmann, Kliimper and Seel have derived novel integral formulas for the correlation 
functions of the spin-1/2 Heisenberg chain at finite temperature. We have found that the high tem- 
perature expansion (HTE) technique can be effectively applied to evaluate these integral formulas. 
Actually, as for the emptiness formation probability P(n) of the isotropic Heisenberg chain, we have 
found a general formula of the HTE for P(n) with arbitrary n £ Z>2 up to 0((J/T) 4 ). If we fix 
a magnetic field to a certain value, we can calculate the HTE to much higher order. For example, 
the order up to 0((J/T) 42 ) has been achieved in the case of P(3) when h = 0. We have compared 
these HTE results with the data by Quantum Monte Carlo simulations. They exhibit excellent 
agreements. 
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The spin-1/2 Heisenberg chain has been one of the 
most fundamental models in the study of the low di- 
mensional magnetism, partially because it can be solved 
exactly by Bethe ansatz. In fact, many physical quan- 
tities of the model have been evaluated exactly even at 
finite temperature [l|. However, they are usually the bulk 
properties, which may be derived directly from the free- 
energy of the system. The exact evaluation of the corre- 
lation functions at finite temperature, on the other hand, 
has remained to be a much more difficult problem. Actu- 
ally the evaluation of the correlation functions is still not 
solved fully even in the case of the ground state, although 
there are various developments recently Note, 
however, that it is established that the correlation func- 
tions in the ground state are expressed in terms of the 
multiple integrals. More recent researches try to evaluate 
these integrals, thereby to obtain theprecise numerical 
values of the correlation functions 0, S B S El ■ 

Considering these situations, it was a great surprise 
that Gohmann, Kliimper and Seel 0, 0, [n| succeeded 
in generalizing the multiple integral formulas to the finite 
temperature, by combining the algebraic Bethe ansatz 
tec hniq ue and the quantum transfer matrix approach 
HI EE 13 13 El- Their results will be a basis for the 
further study of the model at finite temperature. Then 
naturally, it is a significant problem to explore these in- 
tegrals and find a way to extract numerical values of 
the correlation functions at finite temperature. Unfortu- 
nately it is not a straight forward task to generalize those 
methods developed in the case of the ground state to the 
finite temperature. Because the latter formula contains 
an additional auxiliary function a(v), which is a solution 
of a certain non- linear integral equations (NLIE), and is 
more complicated. The purpose of this letter is to at- 
tack this challenging problem with the high temperature 
expansion (HTE) technique. Surprisingly enough, once 
we introduce the HTE, we can perform the multiple in- 
tegrals for each term in the HTE series simply by taking 
a residue at the origin. 



The Hamiltonian of the spin-1/2 isotropic Heisenberg 
chain in a magnetic field h is defined by 
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where a* (k = x,y,z) are the Pauli matrices a k act- 
ing non-trivially on the j-th site of the chain of length 
L. Here we adopt the periodic boundary condition 
(Tj +L — (Tj. In this letter, we mainly consider a spe- 
cial correlation function called the emptiness formation 
probability (EFP) P(n), which is the probability of n 
adjacent spins being aligned upward : 



P(n) = 



Tr{e-f n ;u^} 



Tre 
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At zero temperature T = 0, it was introduced in Q 
and studied further, for example, in 0,0,0. Recently 
Gohmann, Kliimper and Seel obtained [ill Il2j the mul- 
tiple integral formulas of the EFP at finite temperature 



P(n) = 



n 



dyj 



\Jc 2n(l + a( yj )) 
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where functions a(v) and G(v, £) are solutions of the 
NLIE: 
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Here the contour C surrounds the real axis anti-clockwise 
manner. 

First, let us calculate the HTE of a(v) from the NLIE 
Q. This is done by a similar procedure in ^7|> where a 
certain order of the HTE for the free-energy was calcu- 
lated from a NLIE. We assume the following expansion 
for small J/T, 



a(v) = exp 




(6) 



Substituting JJ]J into Q, and comparing coefficients of 
(J/T) m on both sides, we obtain an integral equation 
over {afc(i>)}fcLi for each m (m £ Z>i). As the resultant 
integral equation is linear with respect to a m (v), we can 
solve it recursively. For example, we obtain 



h 2i 
ai(v) = -— - 



a 2 (w) 



J v(l + v 2 ) } 
h 2iv 



J(l + v 2 ) (1 + w 2 ) 2 ' 
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Note that only ai(v) has a pole at the origin. Next let 
us consider the integral equation JSJ. We assume the 
following expansion for small J/T, 



g(«,o=S>(«,o 



fc=0 
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In a similar way, we can determine the coefficients g^{v, £) 
successively by using the results on (0J. For example, we 
obtain 



9o(v,0 = 
9i{v,Q = 



(i+(«-er)(«-o 
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h(2 + 2v 2 - 2vj + £ 2 ) 

'2j{i+v 2 )(i + (v-i) 2 )(i+ey 
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Note that only go(v,£) has a pole at v = £. Finally, 
substituting |JBJ and |(SJ| into we can obtain the HTE 
of P(n). Unexpectedly we have found that we only have 
to calculate residues at the origin. In fact, in this way, 
we could calculate the HTE of the P(n) for small n (n e 
{2,3,4,5,6}). The result up to the order 0((J/T) 4 ) is 



compactly presented as 

1 -2J(-l + n) + hn 



P(n) 



2" 



+ |4J 2 (-4 + n)(-l + ?i) + /i 2 (-l + n)n 
-4M2+(-l + n)n)}— 

+ |l2/iJ 2 (-2 + n) 2 (-l + n) + h 3 (-3 + n)n 2 
-8J 3 (-24 + (-9 + ra)(-3 + n)n) 
-6h 2 J(-l + n)(2 + (-1 + n)n)} 3 _ 2 l_ nT3 

+ |24/ ) , 2 J 2 (-2 + n) 2 (-l + n) 2 

+16J 4 (-5 + n)(-32 + n{26 + (-17 + n)n)) 
+h A n(2 + ra(3 + (-6 + n)n)) 
-32/iJ 3 (24 + (-9 + n)n(6 + (-3 + n)n)) 

-8/i 3 J(-4 + (-4 + n)n(3 + n 2 ))} - ■ 
+ O (( W 5 ) • (10) 
We observe that P(n) has the following form 



2n ^ m! 

m=0 \fe=0 



where po,o( n ) = 1 an< i Pm.fcW for m G Z>i are functions 
of n, which are independent of J, h and T. If we admit 
that Pm.k(n) is a polynomial of n whose degree is at most 
to, our formula l|10(l is also valid for any n S Z>2 as the 
rn-th order polynomial is determined by distinct to + 1 
points. 

Next we fix the magnetic field h to certain values and 
calculate the HTE to much higher order. Then we have 
succeeded in obtaining coefficients of P(3) up to the order 
0((J/T) 42 ) in the case of h = 0. For a finite h, we can 
calculate them at least up to the order 0((J/T) 30 ). It 
will not be easy to obtain HTE coefficients, in particular 
under the magnetic field, to such a high order by other 
method except for the free models. We list some of our 
results on h = case in Table H] 
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FIG. 1: P(n) for J > at h = 0. 
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FIG. 2: P(n) for J > at h = in low T. The exact values 
at T = were obtained in 




FIG. 3: P(n) for J < at ft = 0. 
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FIG. 4: P(3) for J > at h = 0, 2 J, 47, 6 J, ±8J. 
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FIG. 5: P(3) for J < at h = 0, ±0.5] J|, 2| J\, 4| J|. 



Moreover we have applied the Pade approximation to 
our HTE and plot the results in Fig. 1-Fig. 5. Note that 
although the formula (0) was originally derived for J > 
case, our HTE results can be analytically continued to 
J < case. For comparison, we have also performed 
Quantum Monte Carlo simulation (QMC) using recent 
open source softwares in ALPS project [2£|. Especially 
we have chosen the SSE algorithm [2]J so as to treat 
finite magnetic field cases. We have performed the sim- 
ulations with the system size L — 128. In Fig. 1-Fig. 5, 
these QMC data are represented by solid triangles, which 
show excellent agreements with the HTE results. Dis- 
crepancy appears only in the very low temperature re- 
gions, where even the Pade approximation of the HTE 
ceases to converge. We omit these apparent deviations 
of the Pade approximation in Fig. 2. We remark that we 
have also tested the validity of our general formula l|10|) 
up to n — 20. 

In the case of J > and h = 0, we see that -P(n) 
monotonously increases as the temperature increases. 
On the other hand, it decreases monotonously for 
J < 0. In this case we have found P(n) — * l/(n+ 1) 
as T — > 0. Another interesting observation in Fig. 0] 
is that, when J > 0, a peak appears for positive 
values of the magnetic field. Its position moves 
from T = oo to as h increases. For example, 
the peak position T max and the peak P(3) max are 
given by (h/J,T ma */J,P(3) ma *) = (2,12.030,0.13015), 
(4,3.7467,0.18973), (6,1.6904,0.33416), respectively. 
Note that in this case, the critical field is h c = 8 J at 
T = 0, where all the spins are directed upward. 

In conclusion we have found that our HTE method is 
very powerful to evaluate the integral formula for P(n) 
at finite temperature. As an alternative way, it may be 
possible to solve the NLIEs (@J and JHJ) numerically and 
perform numerical integration for the multiple integrals 
in ©. We have tried it, but found it hard to get reliable 
numerical results even for P(3). 

It is straight forward to generalize the results in this 
letter to more general correlation functions. Actually as 
for the nearest and the next-nearest-neighbor correlation 
functions for h = 0, we can immediately calculate their 
HTEs from our results through the relations {SjSj +1 ) = 
P(2)-l/2, (5|5| +2 ) =2(P(3)-P(2) + 1/8), from which 
we can obtain the coefficients whose order is higher than 
the ones by other method 0] . We can evaluate the HTE 
for other complicated correlation functions based on the 
multiple integral formula on the density matrix of the 
XXZ chain at finite temperatures frjl. We will report 
on details in the forthcoming paper |22| . 
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TABLE I: Series coefficients pk(n) for the high temperature expansion of P(n) =J2 k Pk(n)U) k at ft = 0. 
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